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Abstract 

We argue that clustering of color sources, leading to the percolation transition, may 
be the way to achieve deconfmement in heavy ion collisions. The critical density for 
percolation is related to the effective critical temperature of the thermal bath asociated 
to the presence of strong color fields inside the percolating cluster. We find that the 
temperature is rapidity, centrality and energy dependent. We emphasize the similarities 
of percolation of strings with color glass condensate. 
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Recently [1][2] it has been proposed a new thermalization scenario for heavy ion 
collisions which at sufficiently high energies inplies the phase transition to the quark- 
gluon plasma. The key ingredient is the Hawking-Unruh effect [3] [4]. It is well known 
that the black holes evaporate by quantum pair production, and behave as if they have 
an effective temperature of 

T H = k/2tt (f) 

where k = (4GM) -1 is the acceleration of gravity at the surface of a black hole of mass 
M. The rate of pair production in the gravitational background of the black hole can 
be evaluated by considering the tunneling throuh the event horizon. The imaginary 
part of the action for this clasically forbidden process corresponds to the exponent of 
the Boltzmann factor describing the thermal emision [5]. Unruh showed that a similar 
effect arises in a uniformly accelerated frame, where an observer detects an thermal 
radiation with the temperature 

T v = a/27T (2) 

where a is the acceleration. The event horizon in this case emerges through the exis- 
tence of causally disconnected regions of space-time [6] . 

Let us now turn to hadronic interactions. The probability to produce states of 
masses M due to the chromo-electric fiel E and color charge g is given by the Schwinger 
mechanism [7] 

W M ~ exp{-irM 2 /gE) ~ exp(-M/T) (3) 
which is similar to the Boltzmann weight in a heat bath with an effective temperature 

T = |- , a = 2gE/M. (4) 

The full probability P to produce any state of invariant mass M is 

P{M) = 2nW M p(M), (5) 

being p the density of hadronic final states. In the dual resonance model, p is given by 

p{M) ~exp (Wjn , (6) 
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where b is related to the string tension a by the relation b = Xjliso. Unitarity implies 
that the sum of P(M) over all finite states should be finite, what means 



The quantity of the right hand side of (7) is the well known Hagedorn temperature, 
the limiting temperature of hadronic matter above which hadronic matter under goes 
a phase transition into the deconfined phase [8]. So, it also can be considered as 



field. The existence of a limiting temperature is not exclusive of the dual resonance 
model, but it occurs in any system whose exponential mass spectrum has combinatorics 
structure [9]. 

Heavy ion collisions are currently described in terms of color strings stretched be- 
tween the nucleons of the projectile and the target, which decay into new strings 
through g— q pairs production and subsequently hadronize to produce observed hadrons 
[10]. Color strings may be viewed as small discs in the transverse space, 7rrQ,r — 
0.2/m, filled with the color field created by the colliding partons. Particles are pro- 
duced by the Schwinger mechanism, emitting q — q pairs in this field. With growing 
energy and/or atomic number of colliding nuclei, the number of strings grows and they 
start to overlap forming clusters. At a critical density a macros copic cluster appears 
that marks the percolation phase transition [11]. The strong color field inside this 
large cluster produces an acceleration which can be seen as a thermal temperature T, 
by means of the Hawking- Unruh effect. 

In this paper, we compute this temperature to check if its value is above tre critical 
phase temperature found in lattice QCD [12] which is 170-185 MeV depending on the 
number of flavour (2 or 3) and for vanishing chemical potentials. 

In two dimensional percolation theory the relevant variables is the density 




(7) 



the existence of a limiting acceleration a = J 3 /2b corresponding to a strong color 




rj = N s S 1 /S A 



(8) 



where N$ is the number os strings formed in the collisions and S\ and Sa are the 
transverse area of one string an the nuclear overlapping area respectively. Thus, it 
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depends on the impact parameter b of the collision. For very central collisions b = 
and Sa = nR\. 

The critical point for percolating is rjc = 1.18 — 1.5, depending on the type of the 
employed nuclei profile functions, homogeneous or Wood-Saxon [13]. 

We assume that a cluster of N strings behave as a single string with energy momen- 
tum that corresponds to the sum of the energy-momentum of the overlapping strings 
and with a stronger color field resulting from the vectorial sum of the color charges Qi 
of each individual string. The obtained color field covers the area S N of the cluster. As 
the individual string colors are random oriented in color space, the average Qu . Qij 
is zero, such that if S jv = Si then Q 2 N = NQ\ or Q N = -±=NQi with a color reduction 
factor 1/y/N. In general, Sn > Si and the color reduction becomes jf^f [14]. 

In the framework of the Schwinger model extended to color fields, the particle 
density /in is proportional to the color charge Qn and the < p\ >n is proportional to 
the string tension a N ,{i N ~ Q N , and < cr N . On the other hand, from Gauss 

theorem 



Qn ~ S N a N (9) 

we thus obtain 



IS S 
PLn = YTViT^ 1 and < Pt > n= <Pt>i ( 10 ) 

The ratio of the two equations (10) gives rise to the scaling law 

< Pt > S n //j, n =< Pt >i Sx/fj,! (11) 

valid for all projectiles, targets and energies. This scaling law is in reasonable agree- 
ment with data [15]. A similar scaling is found in the Color Glass Condensate (CGC) 
framework [16]. 

Asymtotically, for a random distribution of Ns strings, 

< (12) 

with 

F (V) = J—— (13) 
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being the color reduction factor. It follows that 

H = F(rj)NsfJ,i and < p 2 T >=< p 2 T > 1 /F(r]) (14) 

Note that a N = ai/F(n). 

The tension of the macroscopic cluster fluctuates around its mean value because 
the chromoelectric field is not strictly constant. Assuming a gaussian form for these 
fluctuations [17] we have for the transverse momentum distributions 



dn 



d 2 p± \ < x 2 > 
which gives rise to the thermal distribution 



J dxexp(-x 2 /2 < x 2 >)exp(^-7r^j (15) 



dn / / 2n . 

wr exp {~nTx^\ (16) 

with 

< x 2 >= 7T < p 2 T > l /F(rj) (17) 

Therefore, 



T = y< Pt >ilpF{ri) (18) 

If we identifiy the percolation transition temperature to the Hagedorn temperature, we 
obtain 



which gives for ^< p\ > x ~ 250MeV, and with n ~ 1.18 - 1.5, T c ~ 200 - 250MeV. 

Similar value has been obtained previously in the framework of the CGC [1],[18]. 
In this case T = Qs/27r where Qs is the saturation momentum. As it has been pointed 
out, the role of the scale Qs in the transverse momentum distributions is played by 
< Pt >i / F{rj) in percolation. Similar phenomenology can be obtained in both schemes 
[19]. 

Most of the considerations presented in the framework of the CGC [18] related to a 
rapid thermolization can be applied to our approach. Hydrodynamics models [20] have 
sucessfully described the collective flow of the produced hadrons and the low pt single 
particle spectra. This sucess has led to the conclusion that the created matter behaves 
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like an ideal liquid with almost negligible viscosity [20] [21]. The existing microscopic 
approaches are not able to explain a fast thermalization. For instance, a transport 
model based on independent strings [22] does not present early thermalization although 
the corresponding energy density profiles at proper time r = lfm/c compare well with 
hydrodynamics assumptions for initial energy density distributions. At this point, the 
interaction among string, forming a large cluster with strong color field inside helps to 
reach a fast thermalization. 

In fact, below the percolation threshold we have clusters with a few number of 
strings and different occupied areas what give rise to different temperatures. Above 
the threshold, due to the high density, there are strong interactions among the strings 
forming essentially only one large cluster. The fragmentation of this cluster give rise 
to a thermal distribution of the produced particles. The interaction among these large 
number of strings can be considered as an strong interaction at partonic level which 
give rise to a color composition and the formation of one cluster. 

Similarly to what happens in the CGC framework, we can ask for the charasteristic 
time over which the chromo-electric field changes. In our case, it would be proportional 
to the momentum scale, i.e. < pt >i /F(r)), which implies a lifetime of the cluster, 
r = 1/T ~ 0.7 — lfm for rj > rjc- As the string density increases r decreases. Similar 
conclusions can be obtained, in a more detailed model as it is discussed in ref [23]. 

In the CGC the temperature depends on the rapidity provides that Qs depends on 
y. This dependence may trigger instability of the system and speed up thermalization 
process. In percolation, the temperature (average transverse momentum) depends also 
on the center of mass rapidity y. In fact, the first equation in (14), making ji — > ^ 
and r/i — > n, can be rewritten in the form 

^ = F( V )N s (A,y b )n (20) 

where A = y — y b , ^ ~ exp(y b ) and rj = (^^) N s (A,y b ). Note that, in general, Ns is 
a function of y and y b while n is not. 

For A ~ 0, the rapidity distribution reflects the scaling behaviour of the parton 
structure functions of the valence strings, their number being proportional to the num- 



7 



ber of nucleoli participants, N part = 2N A : 

A~0 N s (0,y b )^N s (0)~Ni(0)N A (21) 

where N$ is the number of strings in pp collisions. 

For A ~ —y b , we are in the domain of sea partons and sea string contribution, their 
number being proportional to the number of nucleon-nucleon collisions, A^ 3 : 

A ~ -y b N s (-y b , y b ) - N s (y b ) ^ N P s (y b )N 4 A /3 (22) 

/,From equations (20),(21)and (22) it follows 

dn I 

A ~ : — ~ N A — yj§ independent of y b (23) 

y n a 

A ~ —y b : ^ ~ N A , increasing with energy (24) 
dy 

Regarding the pt behaviour we have 

A ~ :< p T >— i] 1 ^ 4 — N A 12 '', independent of y b (25) 

A ~ — y b :, < p T >^ r] 1 ^ 4 ~ A^^ 6 increasing with y b (26) 
In this way, the temperature is larger in central rapidity than in the extremes in 

1 /12 

a factor N A , which for central Au — Au collisions is around 1.6. This difference 
increases with energy. This considerable difference of temperature in different rapidity 
slices generater a viscosity, which is very small. Following [18] , the shear viscosity can 
be estimated by 

— =< VT > A ~< p T > —— = (27) 

n nb\ rjp (r/J 

where we use for n, the number of strings per unit of volume, n = L being L the 

longitudinal extension, L ~ lfm. For rj > rjc, we obtain a value around 1 decreasing 

with the density. 
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